Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN MfeK style file vl.4) 



Galaxy groups in the 2dFGRS: the luminous content of the 
groups 

V.R. Eke 1 , Carlos S. Frenk 1 , Carlton M. Baugh 1 , Shaun Cole 1 , Peder Norberg 2 , 
John A. Peacock 3 , Ivan K. Baldry 4 , Joss Bland-Hawthorn 5 , Terry Bridges 5 , Rus- 
sell Cannon 5 , Matthew Colless 6 , Chris Collins 7 , Warrick Couch 8 , Gavin Dalton 9 ' 10 , 
V^- ■ Roberto de Propris 8 , Simon P. Driver 6 , George Efstathiou 11 , Richard S. Ellis 12 , Karl 
^ ■ Glazebrook 4 , Carole A. Jackson 6 , Ofer Lahav 11 , Ian Lewis 9 , Stuart Lumsden 13 , 
■ Steve J. Maddox 14 , Darren Madgwick 15 , Bruce A. Peterson 6 , Will Sutherland 10 , 
*g Keith Taylor 12 (the 2dFGRS Team) 

1 Department of Physics, University of Durham, South Road, Durham DH1 3LE, UK 
^ , 2 ETHZ Institut fur Astronomie, HPF G3.1, ETH Honggerberg, CH-8093 Zurich, Switzerland 
^v^j ' 3 Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK 

' 4 Department of Physics & Astronomy, Johns Hopkins University, Baltimore, MD 21218-2686, USA 

, 5 Anglo- Australian Observatory, P.O. Box 296, Epping, NSW 2121, Australia 

, 6 Research School of Astronomy and Astrophysics, The Australian National University, Canberra, ACT 2611, Australia 



7 Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Birkenhead, L14 1LD, UK 
^Department of Astrophysics, University of New South Wales, Sydney, NSW 2052, Australia 
9 Department of Physics, University of Oxford, Keble Road, Oxford OX1 3RH, UK 
■ 10 Rutherford Appleton Laboratory, Chilton, Didcot, OX11 0QX 
11 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK 
12 Department of Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 
13 Department of Physics, University of Leeds, Woodhouse Lane, Leeds, LS2 9JT, UK 
14 School of Physics & Astronomy, University of Nottingham, Nottingham NG7 2RD, UK 
• 15 Department of Astronomy, University of California, Berkeley, CA 92720, USA 

Or 

6 

' 2 February 2008 



in 



o 
o 



ABSTRACT 

The 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue of - 29000 ob- 
jects is used to study the luminous content of galaxy systems of various sizes. Mock 
galaxy catalogues constructed from cosmological simulations are used to gauge the 
accuracy with which intrinsic group properties can be recovered. It is found that a 
Schechter function is a reasonable fit to the galaxy luminosity functions in groups of 
different mass in the real data, and that the characteristic luminosity L* is larger for 
more massive groups. However, the mock data show that the shape of the recovered 
luminosity function is expected to differ from the true shape, and this must be allowed 
for when interpreting the data. Luminosity function results are presented in both the 
bj and r-p wavebands. The variation of halo mass-to- light ratio, T, with group size is 
studied in both these wavebands. A robust trend of increasing T with increasing group 
luminosity is found in the 2PIGG data. Going from groups with bj luminosities equal 
to 10 10 h~ 2 L to those 100 times more luminous, the typical &j-band mass-to-light 
ratio increases by a factor of 5, whereas the rp-band mass-to-light ratio grows by a 
factor of 3.5. These trends agree well with the predictions of the simulations which 
also predict a minimum in the mass-to-light ratio on a scale roughly corresponding to 
the Local Group. Our data indicate that if such a minimum exists, then it must occur 
at L < 10 10 h~ 2 L , below the range accurately probed by the 2PIGG catalogue. Ac- 
cording to the mock data, the bj mass-to-light ratios of the largest groups are expected 
to be approximately 1.1 times the global value. Assuming that this correction applies 
to the real data, the mean bj luminosity density of the Universe yields an estimate of 
Q TO = 0.26 ± 0.03. Various possible sources of systematic error are considered, with 
the conclusion that these can affect the estimate of Vt m by up to 20%. 

Key words: galaxies: groups - galaxies: haloes - galaxies: clusters: general - large- 
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1 INTRODUCTION 

The distribution of galaxy luminosities, along with their 
mass-to-light ratios and spatial distribution, represent key 
observations that should be explained by theories of galaxy 
formation. In the ACDM model of structure formation, these 
broad brush empirical characterisations of the galaxy pop- 
ulation result from a complicated interplay between the 
lumpy growth and coalescence of dark matter haloes and 
the radiative cooling, star formation and feedback associ- 
ated with baryons within these clumps. Not only are the 
galaxies signposts to locate the underlying overdensities in 
the dark matter distribution, their properties also provide 
a systematic record of the processes that have taken place 
in those haloes and their progenitors. To develop a deeper 
understanding of the impact of these various processes, it 
is helpful to break down the global constraints referred to 
above and consider their variation with halo mass. 

The determination of the luminosity function of galax- 
ies has received much observational effort (e.g. Blanton et al. 
2001, 2003; Cole et al. 2001; Kochanek et al. 2001; Norberg 
et al. 2002), and recently there have been suggestions that 
the distribution of galaxy luminosities varies between rich 
cluster environments and low density regions (e.g. Zablud- 
off & Mulchaey 2000; Christlein 2000; Balogh et al. 2001; 
de Propris et al. 2003). Furthermore, some authors find ev- 
idence that a Schechter function does not well describe the 
galaxy luminosity function in groups and clusters because of 
an excess of bright galaxies (e.g. Smith, Driver & Phillipps 
1997; Trentham & Tully 2002; Christlein & Zabludoff 2003) 
While these empirical results are intriguing, a systematic 
study with a large set of homogenous groups would add wel- 
come weight to these findings. What should one expect to 
find? This question has been addressed using semi-analytical 
prescriptions (White & Frenk 1991; Kauffmann et al. 1993; 
Diaferio et al. 1999; Benson et al. 2003b). These techniques 
provide a physically motivated ab initio procedure for cal- 
culating the properties of galaxy populations, which can be 
compared with the real Universe in order to constrain the 
physical processes and their implementation within the mod- 
els. Both Diaferio et al. and Benson et al. find that the galaxy 
luminosity functions within haloes of different mass are not 
well described by Schechter functions. Instead, they have an 
excess in the abundance at the bright end which arises from 
the different processes that are important for the growth 
of the large central galaxy relative to the other galaxies in 
the haloes (the satellites). In smaller haloes, where a larger 
fraction of the group luminosity is typically locked up in the 
central galaxy, this produces a more prominent deviation 
from a Schechter function at the bright end of the group's 
galaxy luminosity function. Furthermore, the luminosity at 
which this 'central galaxy bump' occurs, increases with halo 
mass. Testing the predictions from these models represents 
an opportunity to learn about aspects of galaxy formation. 

The mass-to-light ratios (T) of groups represent another 
important clue, the light-to-mass ratio essentially reflecting 
the efficiency with which stars are formed within haloes of 
different mass. In large haloes this should be determined 
by the rate of gas cooling, whereas in small haloes other 
factors that cause energy to be injected into the halo gas 
become effective at disrupting the formation of stars. The 
prediction of semi-analytical models is that mass should be 



converted most efficiently into optical light in haloes of mass 
~ 10 12 h,- 1 M (White & Frenk 1991; Kauffmann et al. 1999; 
Benson et al. 2000). One way in which this can be tested is 
to take a theoretically motivated mass function and match 
haloes with galaxy groups of the same abundance as inferred 
from an empirically determined group luminosity function. 
This reveals the variation of the typical T with mass. Qual- 
itatively, Marinoni & Hudson (2002) find a trend similar to 
that predicted by semi-analytical models, namely a mini- 
mum T at intermediate masses of ~ 10 12 /i _1 M o . Obser- 
vationally, there have been many studies that have directly 
measured T for clusters and large groups of galaxies in vari- 
ous wavebands (R - Carlberg et al. 1996, 2001; Tucker et al. 
2000; V - Schaeffer et al. 1993; David, Jones & Forman 1995; 
Cirimele, Nesci & Trevese 1997; Hradecky et al. 2000; B - 
Ramella, Pisani & Geller 1997; Adami et al. 1998; Girardi 
et al. 2000, 2002; Sanderson & Ponman 2003; Tully 2004). 
Some studies suggest that the mass-to-light ratio is larger 
in bigger systems, whereas other studies find no significant 
variation. 

In addition to the direct relevance of the mass-to-light 
ratio and its dependence on halo mass for studies of galaxy 
formation, T is an important quantity to measure since it 
provides one of the traditional ways to estimate the mean 
matter density of the Universe (e.g. Carlberg et al. 1996). 
This requires an assumption about the relation between the 
mass-to-light ratio of the groups or clusters and the mean 
cosmic value. Traditionally, it has been assumed that the 
value of T for the most massive clusters is representative of 
the universal mean, although for a long time this assumption 
was challenged in biased models of galaxy formation (Davis 
ct al. 1985). At any rate, the value of T for clusters, together 
with an appropriate assumption about its universality and 
knowledge of the total galaxy luminosity density allows a 
simple estimate of the mean mass density for the Universe. 

The construction of the 2dFGRS Percolation-Inferred 
Galaxy Group (2PIGG) catalogue (Eke et al. 2004) from the 
2 degree Field Galaxy Redshift Survey (2dFGRS: Colless et 
al. 2001, 2003) facilitates a direct calculation of the galaxy 
content in a large sample of homogenous groups. The pur- 
pose of this paper is to report the results of a decomposition, 
by group size, of the luminosity functions of galaxies within 
groups, and the group mass-to-light ratio. Mock catalogues 
constructed from cosmological simulations are employed as 
a guide to the accuracy with which these quantities can be 
inferred from the 2PIGG sample. This is a vital step in the 
comparison of a model prediction with the observational re- 
sults. 

Section 2 contains a brief description of the group cata- 
logue and a quantitative study of its reliability. In Section 3 
the luminosity functions of galaxies within groups is shown 
for different halo masses. The variation of total halo mass- 
to-light ratio with halo size is investigated in Section 4. 



2 INFERRED GROUP PROPERTIES 

A detailed description of the construction of the mock cata- 
logues, and the groupfinding algorithm applied both to these 
and the real 2dFGRS was given by Eke et al. (2004). Briefly, 
three different ACDM dark matter N-body simulations were 
used. They had (L, N p , <j 8 ) = (154, 288 3 , 0.7), (250, 500 3 , 
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Figure 1. The number of groups at z < 0.12 and z < 0.05 as a 
function of the group dynamical mass, M. Solid and dotted lines 
correspond to the mock catalogues with z max = 0.12 and n m i n = 
2 and 10 respectively. Filled squares and crosses are likewise for 
the 2PIGG catalogue. The dashed line is for z max = 0.05 and 
«min = 2, and the corresponding result from the real catalogue is 
shown with stars. 



0.8) and (141, 256 3 , 0.9), where L is the side of the compu- 
tational cube in h^ 1 Mpc, N p the number of particles and 
as the rms linear density fluctuation in Mpc spheres. 

(The last of these is the GIF simulation described in Jenk- 
ins et al. 1998.) The semi-analytical model of galaxy forma- 
tion described by Cole et al. (2000) was implemented in the 
simulations to create populations of model galaxies whose 
observable properties are given by the model. 

Although the semi-analytical model produces a galaxy 
luminosity function which is quite similar to that in the 
2dFGRS, a small rescaling of the bj luminosities was ap- 
plied, to produce a model luminosity function which was 
identical to that in the 2dFGRS. Mock 2dFGRS catalogues 
were then constructed with the same geometry and position- 
dependent flux limit as the 2dFGRS, and samples of groups 
were selected by applying exactly the same cluster-finding 
algorithm used to generate the 2PIGG catalogue. The tf 
band magnitudes were rescaled in the same way as those in 
the &j band, preserving the semi-analytical galaxy colours. 
Unlike in the bj band, the resulting galaxy rF-band lumi- 
nosity functions for the mock and 2dFGRS catalogues are 
thus not identical, as detailed in Appendix A. 

Here, the focus is on describing the pertinent group 
properties and quantifying the accuracy with which they are 
inferred. The latter relies entirely on the mock catalogues, 
for which the true properties of the haloes hosting the groups 
are known in the parent simulation. 

As described by Eke et al. (2004), the group mass is 
inferred dynamically according to 
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Figure 2. The number of groups as a function of the total group 
bj band luminosity, Lj,j- The lines and symbols have the same 
meaning as in Fig. 1. 
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where A = 5.0, a is the 1-dimensional velocity dispersion, 
calculated using the gapper algorithm (Beers, Flynn & Geb- 
hardt 1990) and removing 85kms _1 in quadrature to ac- 
count for redshift measurement errors, and r is the r.m.s. 
projected separation of galaxies from the group centre, as- 
suming an Q m = 0.3, Qa = 0.7 cosmological model. The 
value of A has been chosen by reference to mock catalogues 
in which galaxies trace the dark matter within each halo, 
apart from a 'central' galaxy placed at the halo centre of 
mass. There is some evidence from gravitational lensing 
studies that bright galaxies do indeed trace the mass, at 
least in rich clusters (Hoekstra et al. 2002; Kneib et al. 2003). 
However, these current observational constraints are not yet 
sufficiently stringent that the effective value of A can be 
empirically justified to better than a few tens of per cent. 
This is the main potential source of systematic error in this 
paper. 



2.1 Calculating group luminosities 

In order to calculate the total group luminosity, it is neces- 
sary to correct for the incompleteness in the 2dFGRS. When 
finding the groups, galaxies from the parent catalogue that 
have no measured redshift have their weights redistributed 
equally to the nearest 10 projected galaxies with measured 
redshifts, and these galaxies are assigned a larger linking vol- 
ume accordingly. As the sets of galaxies with and without 
redshifts are random subsamples (in terms of the intrinsic 
galaxy properties) of the same underlying galaxy popula- 
tion (ignoring the small level of flux-dependent redshift in- 
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Figure 3. The median accuracies (symbols) ±34 percentiles (error bars) of the inferred group properties as a function of the number of 
galaxies in a group, n gal . All groups with z < 0.12 are included. The symbols correspond to the following group properties: y = M (open 
squares), (crosses) and L rp (filled squares). The points have been slightly displaced either side of the true bin values for clarity. 



completeness), the incompleteness can simply be taken into 
account by totting up the observed group bj luminosity ac- 
cording to 



(2-2) 



where the sum extends over the n ga i galaxies in the group 
with their individual weights, u>i. The total group luminosity 
is then estimated by including the contribution from galaxies 
below the luminosity limit, L m i n , at the redshift of the group. 
The extrapolation to zero galaxy luminosity is performed 
assuming that a Schechter function describes the number of 
galaxies as a function of luminosity. Given that the Schechter 
function is 

<KL)dL = 0,(-) exp (--)-, (2.3) 

this operation merely involves dividing the observed lu- 
minosity by the incomplete Gamma function T(a + 
2, L m i n /L»)/r(a + 2). In practice, a small correction is 
also applied because galaxies with fluxes corresponding to 
6j < 14 have been removed from the redshift catalogue. This 
bright flux limit only makes a perceptible difference for a few 



very local groups. These calculations require the additional 
information that M = 5.33 in the bj band and the adopted 
k + e-correction is similar to that of Norberg et al. (2002) , 



k + e 



(2.4) 



l + 8.9z 2 - 5 ' 

The global (M*,a) = (-19.725,-1.18) values are used to 
extrapolate to the total luminosity for all groups. These val- 
ues differ slightly from those of Norberg et al. (2002) as 
they have been derived from the recalibrated 2dFGRS (Col- 
less et al. 2003) . Using a halo mass-dependent extrapolation 
(see Section 3) changes the inferred luminosities by no more 
than ~ 10 per cent. 

At z > 0.12, the fraction of the total group luminosity 
that is actually observed drops below a half. Additionally, 
according to the mock catalogues, the amount of contamina- 
tion of group membership increases at these higher redshifts. 
Thus, in all of what follows, only groups at z < 0.12 will be 
considered. This provides a large number of well-sampled 
groups with a range of masses up to ~ 10 15 h^ 1 M Q . The 
number of groups as a function of mass and luminosity is 
shown in Figs 1 and 2 for both the mock and real catalogues, 
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Figure 4. The median accuracies (symbols) ±34 percentiles of the inferred group properties as a function of the inferred group luminosity, 
L(,j . In the top panel, n m ; n = 2 and, in the bottom panel, ra m i n = 10. Each bundle of three symbols has the following order from left to 
right: y = M, and L rw . The points have been slightly displaced either side of the true bin values for clarity. 



choosing the minimum group membership to be n mln = 2 
or 10. Distributions of n m i n = 2 groups at z < 0.05 are also 
shown in these figures. The mock and real catalogues yield 
broadly similar numbers of groups with similar masses and 
luminosities. There is a slight deficit of groups in the mock 
catalogue relative to the 2PIGG data. This originates from 
the lack of low luminosity galaxies included at z < 0.04, as 
described by Eke et al. (2004). 



Complications arise when using the SuperCOSMOS 
(Hambly et al. 2001) r F band data because the 2dFGRS 
is a bj selected survey. Different depths in the red luminos- 
ity function are thus probed for galaxies of different spec- 
tral types. Rather than dealing with the explicit dependence 
of the luminosity function on colour, the sample was re- 
cut to a conservative r F limit such that almost all 2dF- 
GRS galaxies can be detected to the same red luminosity 
at a given redshift. This issue is discussed in detail in Ap- 
pendix A, which justifies re-cutting the sample at a local 
limit of rF.iim = bj.lim - 1.5. 



2.2 The accuracy of measured halo properties 

The typical redshift error in the 2dFGRS of ~ 85kms _1 
(Colless et al. 2001) will lead to small mass haloes having 
considerable errors in their estimated velocity dispersions 
and hence masses. For instance, a = 85kms _1 and r — 
100 h^ 1 kpc correspond to a mass of about 10 12 h' 1 M Q . In 
order to interpret the results in the following sections, it 
is important to quantify the degree of uncertainty in the 
estimated halo masses and luminosities and how this varies 
with properties such as the number of member galaxies, n ga i, 
or halo luminosity. 

The 'accuracy' with which a halo property is measured 
is defined as a logarithmic measure of the error in the prop- 
erty: 



accuracy = logj 



(^recovered \ 
ytruo / ' 



(2.5) 



where y can stand for mass, &,j or r F luminosity. Figs 3 and 4 
show how the median accuracies, and the 16th and 84th per- 
centiles (as opposed to the mean ±lcr), vary as functions 
of n ga i and Lh 3 respectively. Fig. 3 shows that the spread 
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in mass accuracies becomes very large when the number of 
group members is less than 10. For more populated groups, 
the median accuracy is approximately zero, meaning that 
the typical inferred mass is unbiased. The luminosities are 
typically overestimated by a few tens of per cent even for 
well-populated groups, but they can still be relatively well 
determined for binaries. From this figure, one can infer that 
measured group luminosity is a better quantity than the 
measured dynamical mass for ranking the groups in order of 
size. Note that the accuracy of luminosities has a larger 
spread than that in the bj band. This is the result of using 
fewer galaxies to calculate the luminosity, meaning that 
the Schechter function correction for unobserved galaxies is 
slightly larger in the tf than in the b,j band. Similar trends 
are visible in Fig. 4, where the distributions of accuracies 
are shown as a function of group luminosity for two dif- 
ferent values of n m i n . At low group luminosity, particularly 
for n m i n = 2, the halo masses become increasingly overes- 
timated and poorly determined due to the impact of the 
2dFGRS redshift errors and contamination. Also apparent 
is a tail of underestimated group tf luminosities at low bj 
group luminosity. These are binary groups with no members 
that satisfy the adopted flux limits, and are thus assigned 
zero rF-band luminosity. 

Armed with this quantitative understanding of the 
group catalogue, it is now appropriate to see what infor- 
mation can be extracted concerning the galaxy populations 
in different sizes of halo. 



3 GALAXY LUMINOSITY FUNCTIONS 
WITHIN GROUPS 

This Section addresses firstly the issue of how well the galaxy 
luminosity function in haloes of different mass can be recov- 
ered from the mock catalogues and, secondly, the determi- 
nation of these functions for the real data. In calculating 
the space density of galaxies, a 1/V ma x estimator has been 
applied to each galaxy. The variable flux limit across the 
survey has been taken into account. Dynamically inferred 
group masses, calculated according to equation (2.1), have 
been used to split the total sample into different classes. 



3.1 bj band results 

Sets of mock catalogues have been constructed from two dif- 
ferent semi-analytical models of galaxy formation applied to 
the same N-body simulations. Both these models are based 
on the general scheme developed by Cole et al. (2000), but 
they treat certain physical processes in different ways. The 
first model, described by Benson et al. (2002) , is very simi- 
lar to the original Cole et al. model except that it includes 
detailed treatments of photoionization of the intergalactic 
medium at high redshift and of the dynamics of satellites in 
haloes. This model, which will be referred to as the 'bumpy' 
model for reasons that will become apparent shortly, was 
used by Benson et al. (2003b) to predict the luminosity 
function of galaxies in haloes of different mass. The sec- 
ond model comes from Benson et al. (2003a) and will be 
referred to as the 'superwind' model. This has a baryon 
fraction twice as high as that assumed by Cole et al., in 



accordance with recent determinations, and includes a treat- 
ment of superwinds. These suppress the growth of the overly 
bright galaxies which otherwise tend to form in high baryon 
fraction models. 

Fig. 5 shows galaxy luminosity functions in groups of 
different mass in both the semi-analytical models of galaxy 
formation. The solid lines correspond to the bumpy model 
while the dashed lines correspond to the superwind model. 
The suppression of bright central galaxies in the superwind 
model is apparent in the middle and right-hand panels which 
show the galaxy luminosity functions within groups of mass 
M ~ 10 14 /i _1 M o and M ~ 10 15 h' 1 M respectively. It 
is clear that the detailed shapes of the galaxy luminosity 
functions in groups depend upon the assumptions that go 
into making a semi-analytical model. These data can thus 
be used to test the models. In what follows, the bumpy mock 
catalogues will be shown unless otherwise stated. 

In order to recover faithfully both the amplitude and 
shape of the galaxy luminosity function in groups of dif- 
ferent mass, it is helpful to adopt separate values of n m in 
and z max for the various group samples. For instance, 
clusters that should contribute to the highest mass bin 
will have large numbers of members and be detectable 
to higher redshifts compared to the smaller groups. Thus, 
the following empirical choices of (n m i n , z m ax) were made 
for the 10 13 , 10 14 and 10 15 h' 1 M Q samples respectively: 
[(3,0.08), (10,0.12), (80,0.12)]. These will be used through- 
out this Section. 

The results of these selections can be seen in Fig. 6, 
which shows the semi-analytical model luminosity functions 
in the simulations from which the mock catalogues were 
created, the luminosity functions actually recovered from 
the mocks, and the luminosity functions estimated from 
the 2PIGG sample for the same values of n m in and z max . 
While the theoretical luminosity functions in the simulations 
from which the mock catalogues were made differ substan- 
tially from Schechter functions, the bumps have been largely 
smeared out in the mock recovered luminosity functions, 
which now better resemble Schechter functions. This smear- 
ing is predominantly the result of contaminating galaxies 
contributing to haloes with inappropriate masses and ho- 
mogenizing the samples. Note that the original bumps oc- 
cur at luminosities that vary with halo mass. There is a 
tentative bump detection in the two higher mass group sam- 
ples. This excess of luminous galaxies is sufficient to render 
the Schechter function a bad fit, although it does provide 
a good description of the results for the M ~ 10 13 h^ 1 M Q 
groups. Overall however, apart from the smearing out of 
the central galaxy bumps, the shapes and amplitudes of the 
group galaxy luminosity functions are recovered well from 
the mock catalogue. 

Having determined the efficiency with which group 
mass-dependent galaxy luminosity functions are recovered, 
it is now appropriate to consider the real data. The 2PIGG 
results are shown by the points in Fig. 6. Error bars are the 
Poisson errors on the number of groups contributing galax- 
ies to each luminosity bin. In almost all respects, the real 
2PIGG results look very similar to those recovered from the 
mock catalogue, shown as dotted lines. The main difference 
is the abundance of low luminosity (L < L*/3) galaxies in 
higher mass haloes. The mock groups contain more of these 
galaxies than are present in high mass 2PIGGs. This differ- 
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Figure 5. Galaxy 6,j luminosity functions in haloes of different mass predicted by two semi-analytical models of galaxy formation. The 
panels refer to haloes with masses of 10 x h- 1 M Q , with x = 13 ±0.5 (left), 14 ±0.5 (centre) and 15 ±0.5 (right). The scales at the top of 
the panels give the absolute fej magnitude. Solid lines correspond to the 'bumpy' semi-analytical model, whereas dashed lines correspond 
to the 'superwind' model. 
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Figure 6. Galaxy 6j luminosity functions in groups of different mass in both mock and real 2PIGG catalogues. The group mass ranges are 
M ~ 10 13 h^ 1 M (left), M ~ 10 14 fr 1 M Q (centre) and M ~ 10 15 ft. -1 Mq (right). In each case, a solid line shows the semi-analytical 
prediction from which the mock catalogues were constructed, a dotted line shows what is actually measured in the mock groups, and 
the filled squares with errors depict the results from the 2PIGG data. The values of n m i n and 2 m ax are selected for each mass range, as 
described in the text, in order to optimize the recovery of the original galaxy semi-analytical luminosity functions. 



ence is apparent in the panels showing the M ~ 10 14 bT 1 Mq 
and M ~ 10 15 K~ x Mq results. The evidence for bumps in 
the 2PIGG data is also only tentative. 

Fig. 7 shows the 1, 2 and 3 — a contours in the L* — a 
plane for the STY (Sandage, Tammann & Yahil 1979) fits 
to the 2PIGG galaxy luminosity functions within groups of 
different mass. Also shown, without a symbol marking the 
most likely value, are the contours representing the best- 
fitting parameters for the whole galaxy population out to 
z = 0.12. The parameters for the highest mass bin are much 



less well constrained because there are fewer total galax- 
ies in these groups than in the more abundant, lower mass 
haloes. The 2PIGG sample has a slightly brighter L» for 
the M ~ 10 14 h' 1 M groups than the M ~ 10 13 h' 1 M 
groups. This trend is reproduced in the bumpy mock cata- 
logue but not in the superwind mock, in which the most lu- 
minous galaxies no longer reside in the most massive haloes. 

Fig. 8 shows how the recovered luminosity functions dif- 
fer from the suitably normalised best-fitting Schechter func- 
tions for the 2PIGGs. While the fit works well for both mock 
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Figure 8. The ratios of galaxy luminosity functions to the 2PIGG best-fitting Schechter functions for the M ~ 10 13 h^ 1 Mq (left), 
M ~ 10 14 ft, -1 Mq (centre) and M ~ 10 15 ft, -1 Mq (right) groups. Results are shown for 2PIGG (squares), the bumpy mock (dotted 
lines) and the superwind mock (dashed lines). 
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Figure 7. 1, 2 and 3 — <r contours showing the relative probability 
that different Schechter function parameters, L„ and o, provide a 
good description of the ft j-band galaxy luminosity function within 
2PIGGs of different mass. These results were obtained by the STY 
estimation method using the groups contributing to the 2PIGG 
luminosity functions in Fig. 6. The different group masses have 
different symbols marking the most likely parameter values, as 
detailed in the figure. The small ellipses with no central symbol 
represent the results for all 2dFGRS galaxies at z < 0.12. 



is found in the 2PIGGs. Finally, the lack of very luminous 
galaxies in the superwind mock shows that, despite the im- 
perfections of the recovery, it is still possible to discriminate 
between this and the bumpy model. 



3.2 tf band results 

The same analysis performed in Section 3.1 for the bj band 
data can be performed for the band data. While all galax- 
ies can be used to calculate the luminosity functions using 
the 1/Knax method, the STY estimation of the best-fitting 
Schechter function parameters requires a complete sample, 
so galaxies are only included if 14 < < &j,n m — 1.5 when 
determining L* and a. 

Figs 9, 10, 11 and 12 are the tf band equivalents of the 
previous four figures for the bj band. Apart from larger er- 
ror ellipses on the recovered Schechter function fits, resulting 
from the smaller number of galaxies being used to ensure a 
complete sample, the rF band results are qualitatively very 
similar to those found in the bj band. The 2PIGGs with 
again have an average galaxy luminos- 
ity function with L* about 25 per cent higher than that for 
the M ~ 10 13 h' 1 Mq groups, and a slightly steeper faint 
end slope. The bumpy mock reproduces this scaling of L» 
with halo mass, but again contains an excess of low lumi- 
nosity galaxies in the most massive haloes. Once more, the 
superwind mock exhibits both an excess of low luminosity 
galaxies and a deficit of high luminosity galaxies in both of 
the higher mass bins, relative to the 2PIGG sample. 



and real data in the M ~ 10 13 h~ 1 Mq groups, it becomes 
an increasingly poor description in more massive systems, 
where an excess of luminous galaxies can be seen to distort 
the fit. Also apparent is the overabundance of low luminosity 
galaxies in the most massive mock groups relative to what 



4 GROUP MASS-TO-LIGHT RATIOS 

As outlined in the Introduction, the mass-to-light ratio of 
groups, T, contains clues to the nature of galaxy formation, 
and it can also be used to estimate the mean mass density 
of the Universe. This Section contains a description of the 
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Figure 9. The rp band equivalent of Fig. 5, showing the galaxy luminosity functions in groups of different mass in the two semi-analytical 
models of galaxy formation. As before, the panels correspond to groups with M = 10 13±0 - 5 h- 1 M© (left), M = 10 14±0 - 5 h' 1 M (centre) 
and M = 10 15±0 - 5 ft, -1 Mq (right). Solid lines correspond to the 'bumpy' semi- analytical model, whereas dashed lines correspond to the 
'superwind' model. 
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Figure 10. The rp band equivalent of Fig. 6, showing both how well the galaxy luminosity functions can be recovered in different sized 
groups, and how the mock results compare with the real 2PIGG data. As before, the group mass ranges arc M ~ 10 13 h^ 1 Mq (left), 
M ~ 10 14 hr 1 Mq (centre) and M ~ 10 15 h~ x Mq (right). In each case, a solid line shows the semi-analytical prediction from which the 
mock catalogues were constructed, a dotted line shows what is actually measured in the mock groups, and the filled squares with errors 
depict the results from the 2PIGG data. 



accuracy with which T can be determined from the 2PIGG 
catalogue, and addresses the dependence of T on group lu- 
minosity, the size of halo in which stars form most efficiently, 
and the estimate of the mean mass density, Q m - 

In trying to recover the typical group mass-to-light ra- 
tio for groups of a particular size, there is no reason to 
use the entire group sample. One merely requires an un- 
biased subset of groups. Experimentation with the mock 
catalogues has yielded an appropriate subset, in which the 
bias of the results is minimized. In this scheme, the smaller 



groups are only used when they are relatively nearby and 
isolated, whereas the restrictions are less stringent for larger 
groups. Specifically, of the n m i n = 2 groups, only those with 
z < z max = 0.07 + 0. 02[log 10 (L 6j //i _2 L Q )-10] and having 
no neighbouring groups with centres at distances less than 
rfmin/ h^ 1 Mpc = 2 + [ 10 — log 10 (£[,,/ h~ 2 Lq)] were used to 
calculate the mass-to-light ratios. 

The reasons for the success of this empirical choice are 
as follows. Many of the small groups are fragments cleaved 
from much bigger groups by the group-finding algorithm. 
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Figure 12. The rp band equivalent of Fig. 8, showing the ratios of galaxy luminosity functions in the 2PIGG (squares), bumpy mock 
groups (dotted lines) and superwind mock groups (dashed lines) to the best-fitting Schechter function for the 2PIGGs in each mass 
range. The mass ranges are: M ~ 10 13 h,- 1 M (left), M ~ 10 14 h,- 1 M Q (centre) and M ~ 10 15 h' 1 M Q (right). 
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Figure 11. The rp band equivalent of Fig. 7, showing the 
STY-determined most likely Schechter function parameters for 
2PIGG galaxy luminosity functions. The different group masses 
are represented by the contours enclosing the open square (M = 
10 13±0.5 h -l M0 ) ; cross ( M = 10 14±0.5 h -\ Mq -) and fiUed 

square (M = 10 15±0,5 ft -1 Mq). Results for all galaxies satis- 
fying 14 < rp < fej^im — 1.5 and z < 0.12 are shown by contours 
with no central symbol. 

These typically have large velocity dispersions, reflecting the 
size of the halo to which they really belong. Consequently 
their mass-to-light ratios are unrepresentatively high. The 
nearest neighbour distance restriction is intended to elimi- 
nate these spurious groups from the analysis. The reason for 
the redshift limit is the desire that most underlying groups 



of a particular size should contain at least two galaxies that 
will be detected in a flux-limited survey. The mock cata- 
logues used here are such that the true mass-to- light ratio of 
low luminosity groups is typically lower for groups contain- 
ing only one detectable galaxy. Thus, in order to ensure that 
the low luminosity groups are representative of the under- 
lying distribution, only the nearby low luminosity examples 
are included. 

The following two subsections describe the bj and rp 
band results. The redder band traces the stellar mass more 
faithfully than the blue band, and T rp thus provides addi- 
tional useful information. 

4.1 6,j band results 

Fig. 13 shows the dependence of on group mass. The 
horizontal line indicates the mean mass-to-light ratio in the 
mock universe, and the other solid line represents the vari- 
ation of the median group with halo mass in the sim- 
ulation with semi-analytical model galaxies from which the 
mock catalogues were constructed. A dotted line traces what 
is actually recovered from the mock catalogue, and the filled 
squares show the results from the real 2PIGG subsample of 
groups. Error bars represent the 16th and 84th percentiles 
divided by the square root of the number of groups in each 
bin. As a result of the large spread of accuracies in the in- 
ferred group masses for groups containing only a few mem- 
bers, there is a strong smearing effect along a line of gradient 
one. For the highest mass groups, for which more galaxies 
yield better mass estimates, this trend is reduced, but there 
is sufficient contamination from low luminosity groups with 
greatly overestimated masses, that the recovered T;, T is still 
biased high by 30 — 40 per cent. The main reason why the 
mock recovered variation resembles the true behaviour in 
the parent simulation at intermediate masses is that this is 
where the majority of the groups are found (see Fig. 1), so 
there are comparable numbers of groups biasing the results 



© 0000 RAS, MNRAS 000, 000-000 



Luminous content of 2dFGRS galaxy groups 11 



, 2.5 



1.5 



1 ' ' ' ' I 1 

- Model 

- Mock recovered 
i 2PIGG 



=0.30 




gradient 1 



12 13 14 

log lo [M/(h-'M )] 



15 



Figure 13. The group mass-to-light ratios as a function of group 
mass. The solid line traces the variation of the median Y;, T in 
the semi-analytical model from which the mock group catalogues 
were constructed. The dotted line shows the Tj, actually recov- 
ered from the mock catalogue using groups with n m i n = 2 that lie 
within a redshift of 2 max = 0.07 + 0.02 [log 10 (L bl / h~' 2 L©) - 10] 
and have no neighbouring groups with centres at distances less 
than d min / Mpc = 2 + [ 10 - log 10 (L[ )J / h~ 2 L Q )] . The cor- 
responding measurement from the 2PIGG sample is shown by 
the filled squares, with error bars representing the 16th and 84th 
percentiles divided by the square root of the number of groups 
contributing to each bin. This would be the standard deviation on 
the median if the Tj, values were distributed Normally in each 
bin. The horizontal line indicates the mean mass-to-light ratio 
of the mock universe, calculated using the luminosity function of 
Norberg et al. (2002), and the appropriate value of Q m = 0.3. A 
line with gradient 1 is also plotted to illustrate the direction in 
which errors in the inferred group masses would move the esti- 
mated values. 



high and low. In summary, this figure demonstrates that it is 
unhelpful to plot correlated variables when they both suffer 
from the same large uncertainty. 

A much better method for ordering groups in terms of 
size is to use their measured luminosities; these have much 
smaller errors than the group masses, particularly for groups 
with only two members (see Fig. 3). Now, instead of mass 
errors smearing out the results along a line of gradient one, 
there are luminosity uncertainties that act along a line of 
gradient minus one. The sign of this gradient is opposite to 
the trend of Yf,, with halo luminosity in the parent simula- 
tion, and the amplitude of this error vector is now substan- 
tially reduced. Fig. 14 shows the corresponding variation of 
group mass-to-light ratio with group luminosity. 

In the semi-analytical model, the variation of T bj re- 
flects the efficiency of star formation in haloes of differ- 
ent mass. Star formation is most efficient in haloes with 
L ~ 6 x lO 9 /i~ 2 L0, the point at which the minimum 
of the solid line occurs. The dotted line traces the varia- 
tion of Tbj recovered from the group mock catalogue. It 
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Figure 14. The median mass-to-light ratio of groups as a func- 
tion of group luminosity. Line types and symbols are the same as 
in the previous figure, as are the values of n m ; n , z max and d m j n 
defining the recovered groups. 



is immediately apparent that this recovery, for haloes with 
L&j > 10 10 h~ 2 Lq, is very much better than was the case 
when group mass was used to quantify halo size. In par- 
ticular, both the rise in Yb T from small groups to clusters, 
and the plateau for the largest clusters are well reproduced. 
There is a slight bias low, by ~ 40 per cent, during the 
rise. This can be understood in terms of the position of the 
peak in the distribution of groups as a function of lumi- 
nosity (Fig. 2), and the typical errors being made, whereby 
luminosities are slightly overestimated as a result of contam- 
ination (Fig. 4). The recovery of the mass-to-light ratio for 
groups less luminous than L* (for galaxies) is poor. Among 
the reasons for this is the fact that the individual redshift 
measurement errors become comparable with the typical 
group velocity dispersions in such small systems. Also, there 
is insufficient volume to achieve small statistical errors while 
ensuring that these groups contain more than one detectable 
galaxy. 

For the 111 groups with \og w [L b J{ h~ 2 Lq)] > 11.5 in 
the mock catalogue, the recovered median and mean mass- 
to- light ratios are 471 and 549 ± 28 h Mq /Lq respectively. 
The corresponding values in the parent simulation are 507 
(median) and 520 (mean). Thus, the mean is overestimated 
in the mocks by ~ 6 per cent, while the recovered median 
is 7 per cent too low. As the mean is more easily affected 
by outlying data points with high measurement errors, in 
what follows only the median cluster mass-to-light ratio will 
be used. The estimate of Q m from the cluster mass-to-light 
ratio is based on comparing this to the mean mass-to-light 
ratio of the universe. In the parent simulation, the mean 
universal value is 458 h Mq /Lq , 10 per cent lower than the 
median recovered mass-to-light ratio of clusters. Thus, for 
this mock catalogue at least, the median mass-to-light ratio 
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Figure 15. The median mass-to-light ratio of groups as a func- 
tion of group luminosity. As in the previous figure, filled squares 
represent results for the 2PIGG sample. The horizontal line in- 
dicates the mean mass-to-light ratio determined from the mean 
luminosity density of Norbcrg ct al. (2001), assuming Q m = 0.3. 
The results of four other studies of B-band group mass-to-light 
ratios are also shown, as indicated in the legend. Vertical error 
bars on the data points give the uncertainty on the median Tf, r 
The horizontal error bars represent the 16 — 84 percentile range 
of the group luminosities contributing to these medians. 



of the recovered clusters allows a pretty accurate estimate 
of the mean mass density of the simulated universe. 

The filled squares with error bars in Fig 14 represent 
the real 2PIGG mass-to-light results. These are remark- 
ably similar to what was recovered from the mock cata- 
logue, showing the same three basic features: nothing use- 
ful at L £ 10 10 h- 2 Lq, a rise by a factor of ~ 5 up to 
L ~ 2 x 10 11 /i~ 2 Lq, and a constant mass-to-light ratio 
for the largest groups. Unlike the case when mass was used 
to quantify the group size, this trend of increasing median 
mass-to-light ratio with larger groups can no longer be con- 
fused with the effect of mass measurement errors. For the 
log 10 [Lb J /( h~ 2 Lq)] > 11.5 groups, of which there are 96, 
the median Tf,, = 429 ±25 TiMq/Lq, where the uncertainty 
represents the statistical error on the median. This number 
is used later to estimate Cl m . 



4.2 Comparison with other studies 

When comparing these results with other studies, it is essen- 
tial to bear in mind that there exists a variety of definitions 
of a group. To reiterate, the 2PIGG analysis is designed to 
recover, as well as possible, groups that resemble those iden- 
tified in a ACDM dark matter simulation using the standard 
linking length of 0.2 times the mean interparticle separation. 
Note that in Fig. 15, the raw 2PIGG results are shown for 
the purpose of comparison with other studies, despite the 
fact that they probably suffer from a small bias of the sort 



shown in Fig. 14. In Section 4.4, this bias is corrected and 
functions are given which describe the variation of the cor- 
rected median group mass-to-light ratio with group lumi- 
nosity. 

The results of Girardi et al. (2002; G2002), which ex- 
tend the database used by Girardi et al. (2000) with groups 
mostly found in the Nearby Optical Galaxy (NOG) cata- 
logue, contain groups from a number of different surveys. 
Note that these papers both contain a typo relating to the 
conversion between B and bj magnitudes (Girardi, private 
communication). The correct relation, for a mean galaxy 
colour of (B-V) = 0.9, is 



B-bj = 0.252. 
In conjunction with 
(B-6j)© =0.15, 
this leads to 



(*). 



1.1, 



(4.1) 



(4.2) 



(4.3) 



as was correctly reported in G2002. Fig. 15 shows the data 
from the (CL+PS) sample of G2002, supplemented with 500 
other systems from their set of groups found by applying a 
percolation algorithm to the NOG survey (Girardi, private 
communication). Relative to the homogenous 2PIGG sam- 
ple, this set of groups comes from a number of different 
sources and, furthermore, they were analyzed using differ- 
ent types of mass and luminosity estimators. Yet, a broadly 
similar behaviour of the mass-to-light ratio is seen. Care 
is needed, however, because the percolation algorithm of 
G2002 will presumably yield results that are biased to low 
mass-to-light ratios, for the same reason that the 2PIGG 
results are affected. Their percolation algorithm (Giuricin 
ct al. 2000) employs a scaling of the linking volume with 
redshift that is different from that used for the 2PIGG cat- 
alogue (Eke et al. 2004), so the properties of the recovered 
groups will be different. 

Sanderson & Ponman (2003) used a sample of 32 sys- 
tems with masses inferred from their X-ray emission. These 
groups are concentrated towards the large mass end of the 
distribution and show little or no trend of mass-to- light ratio 
with group luminosity, comparable with the 2PIGG results. 
Again, their results yield slightly lower values of the mass- 
to-light ratio than found for 2PIGG. Note that the point 
plotted in the figure differs very slightly from that in the 
Sanderson & Ponman paper as a result of the correction of 
an error in their B to bj conversion for 3 groups. 

Ramella, Pisani & Geller (1997) used a sample of 406 
groups identified in the northern CfA2 survey, with masses 
estimated from optical data. While the median mass-to-light 
ratio is lower than that of the Sanderson & Ponman groups, 
the CfA2 groups are typically smaller and the 2PIGG re- 
sults indicate that Yj, T should be smaller for such groups. 
However, once more, the mass-to-light ratio in this sample 
is lower than that inferred from the 2PIGG catalogue. 

Marinoni & Hudson (2002) have indirectly inferred a 
mass-to-light ratio variation for groups by finding the map- 
ping between the Press-Schechter (1974) mass function, for 
some assumed cosmological model, and the luminosity func- 
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tion of groups measured from the NOG sample by Marinoni 
et al. (2002). Their results are given in the B band, and 
have been converted to the bj band using equation (4.3). 
This curve shows a similar trend to the 2PIGG results but, 
again, with an offset. In general, estimators of the luminosity 
function of groups tend to be biased towards high luminosi- 
ties as a result of the inevitable inclusion of interlopers in 
the groups, the statistical errors in the inferred group lumi- 
nosities and the steep decline in the abundance of objects 
with increasing size. This would lead to an overestimate of 
the luminosity at a particular abundance (or mass), biasing 
low the inferred mass-to-light ratio. Indeed, estimating the 
luminosity function of groups in mock 2PIGG catalogues 
suggests that this bias could plausibly yield an overestimate 
of the group luminosity by a factor of at least 2 at a par- 
ticular abundance (Eke et al. in preparation). If a similar 
size of bias were present in the sample used by Marinoni & 
Hudson, then this would account for their significantly lower 
mass-to- light ratios at L ~ 10 11 h~ 2 Lq. 

One factor that would bias the 2PIGG mass-to-light val- 
ues slightly high is the incompleteness in the parent 2dFGRS 
catalogue (Pimbblet et al. 2001; Norberg et al. 2002; Cross 
et al. 2004). Such an effect, however, is small compared to 
the difference between the 2PIGG results and those of the 
other studies shown in Fig. 15. The most likely other reason 
why the 2PIGG results could be systematically incorrect is 
if the mass estimator employed here is not appropriate in 
the real world, as discussed in Section 2. This would happen 
if galaxies traced real dark matter haloes in a different way 
from what was assumed in making the mock catalogues, ei- 
ther because they had a different velocity dispersion to the 
underlying dark matter, or because they had a different pro- 
jected spatial distribution (De Lucia et al. 2004; Diemand, 
Moore & Stadel 2004). In order for biases of this kind to 
give rise to the observed variation of the mass-to-light ratio 
with group size, however, they would need to depend very 
strongly on group size. It is difficult to imagine that such 
biases could give rise to the factor of ~ 5 variation seen in 
the data. 



4.3 tf band results 

As described in Appendix A, the 2dFGRS combined with 
the SuperCOSMOS r-p data provides a well defined sample 
of galaxies with 14 < tf < 6j,ii m — 1.5. While applying 
these flux limits reduces the number of galaxies available 
to calculate the group luminosity, the missing faint galaxies 
contribute little to the total, so the accuracy with which 
group tf band luminosities can be inferred is comparable to 
that for the bj band, as was shown in Section 2.2. 

Fig. 16 shows how well the variation of T rp in the par- 
ent simulation can be recovered in the mock catalogue, as 
a function of group bj luminosity. (This is a slightly better 
determined measure of 2PIGG size than the tf luminosity.) 
The results are remarkably similar to those shown in Fig. 14 
for Yf,,. Firstly, the intrinsic variation present in the par- 
ent simulation is quite accurately recovered in the mocks. 
Secondly, as was the case for T bj , the 2PIGG data agree 
remarkably well with the behaviour seen in the mock cata- 
logue. The main difference between the tf and bj results is 
the size of the change in typical mass-to-light ratio over the 
reliably probed range of group luminosities. This is not as 
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Figure 16. The variation of the r band group mass-to-light ratio 
with group 6 j luminosity. A solid line traces the median behaviour 
in the parent simulation, and the dotted line shows what is actu- 
ally recovered in the mock catalogue. Filled squares represent the 
results from the 2PIGG catalogue, using the same redshift and 
minimum group separation as were applied to the 6 j band data. 
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Figure 17. The variation of the rp band group mass-to-light 
ratio with group luminosity. Filled squares represent the results 
from the 2PIGG catalogue, and the other symbols show other 
observational determinations, as detailed in the legend. Note that 
the Carlberg et al. results have luminosities measured in the Gunn 
r band, and the Tucker et al. data are in the LCRS R band. 
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Figure 18. The corrected variation of the 2PIGG group mass- 
to-light ratio with group luminosity, in both the bj and rp bands. 
Equations (4.4) for and (4.5) for T rp are shown as solid 

and dashed lines. The corresponding 2PIGG data are shown with 
open squares and crosses. 



extreme in the tf band, a factor of ~ 3.5 rather than the 
factor of ~ 5 seen in the bj band - as one would expect if 
halo size has a stronger effect on recent star formation than 
on the overall stellar mass. 

Fig. 17 compares the 2PIGG results with those from 
CNOC groups and clusters (Carlberg et al. 1996, 2001) and 
the LCRS groups (Tucker et al. 2000). Unlike in the bj band 
case, there is now some agreement between the 2PIGGs and 
some other work. Note that the Carlberg et al. data points 
have luminosities measured in the Gunn r band and include 
a correction to redshift zero, whereas the Tucker et al. data 
are measured in the LCRS 7? band. Once again, the hori- 
zontal error bars illustrate the 16 — 84th percentile range in 
group luminosity contributing to each point. 



4.4 Corrected mass-to-light results 

It was shown in Sections 4.1 and 4.3 that the values of T 
as a function of L bj recovered in the mock catalogues are 
slightly biased low relative to the values in the parent simu- 
lation. This bias is reproducible in mocks constructed from 
other simulations, for example, from the less clustered mass 
distribution of a ACDM simulation with as ~ 0.7 analyzed 
at z « 0.1 (rather than the default simulation with cr 8 = 0.9 
and z — 0.1). 

It seems plausible that a bias of this kind is also present 
in the 2PIGG sample results. A correction factor is defined 
as the ratio between the luminosity-dependent values of T 
measured in the parent simulation and in the mocks (in both 
the 6j and tf bands). The 2PIGG results are simply multi- 
plied by this factor to obtain an estimate of the true underly- 
ing group mass-to-light ratio variation in the real Universe. 
The corrected values are shown with open squares (bj band) 



and crosses (rF band) in Fig. 18. The error bars on the data 
points still represent the size of the statistical uncertainty 
on the median. Approximate fits to the corrected 2PIGG 
median mass-to-light ratios are shown by the solid (bj) and 
dashed (r^) lines, which have the following functional forms: 



log 10 T bj = 2.28 

+ 0.4tanh{1.9[log 10 (L 6j ) - 10.6]} 

and 

logioTr F = 2.20 

+ 0.27tanh{2.4 [log 10 (L bl ) - 10.75]}. 



(4.4) 



(4.5) 



These equations are valid for 10 < log 10 (L 6j / h~ 2 Lq) < 
12.5. Since the 2PIGG results are so similar to those recov- 
ered from the mock catalogue, these fits to the underlying 
'truth' in the real Universe also provide a decent description 
of the behaviour in the parent simulation, as traced by the 
solid lines in Figs 14 and 16. The correction is sufficiently 
small that all the basic features of the results discussed above 
remain. 

The correction changes the median mass-to-light ratio 
of the 98 clusters with log 10 [Li, I /( h~ 2 L )] > 11.5 from 
T 6j = 427 ± 24/iM o /L to 466 ± 26feM /L Q . Together 
with the global mean fej band luminosity density and its un- 
certainty, inferred from the luminosity function measured by 
Norberg et al. (2002), this implies* that Q. m = 0.31±0.03. If, 
as was the case in the mocks, the median mass-to-light ratio 
of the largest clusters overestimates the mean mass-to-light 
ratio of the Universe by ~ 11 per cent, then this estimate 
should come down accordingly to Q, m — 0.28 ± 0.03. This 
uncertainty only includes the statistical error. 

Possible sources of systematic uncertainty, in addition 
to the two model-dependent corrections already applied to 
the raw 2PIGG results, would come from systematic errors 
in the estimation of the cluster mass-to-light ratios. These 
could arise from the underestimate of the cluster luminosi- 
ties resulting from incompleteness in the parent 2dFGRS 
catalogue. In fact, Norberg et al. (2002) do assume a 9 per 
cent incompleteness in the 2dFGRS parent catalogue when 
computing the mean fej-band luminosity density of the Uni- 
verse, so for consistency the 2PIGG luminosities should be 
similarly increased. The resulting lower mass-to-light ratios 
yield Q, m — 0.26 ± 0.03. Another, less readily quantified, 
potential systematic uncertainty would come from a bias in 
the mass estimator if real galaxies do not trace dark matter 
haloes in the way assumed in the mocks. As discussed in 
Section 2 this could plausible be a few tens of per cent. 

It is interesting that both this estimate of fi m and its 
statistical uncertainty are very similar to the values quoted 
by Spergel et al. (2003) from the combination of WMAP 
microwave background data and the 2dFGRS galaxy power 
spectrum. This agreement, in itself, suggests that the sys- 
tematic uncertainties in either estimate are likely to be 
small and provides a welcome consistency check of the en- 
tire paradigm of structure formation by hierachical cluster- 
ing from CDM initial conditions. Further reassurance can 
be gained from using the rF-band mass-to- light ratios of the 



* Assuming that the galaxies trace the mass in clusters, as is 
customary when using this method 
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same set of the most luminous clusters and the rF-band lu- 
minosity density of the Universe (see Appendix A). This 
yields an estimate of fi m = 0.25 ± 0.04, which is consistent 
with the value inferred from the fej-band data. 



5 CONCLUSIONS 

This paper presents an analysis of the galactic content of 
the 2PIGG catalogue of groups and clusters identified in 
the 2dFGRS, focusing on the galaxy luminosity function and 
total mass-to-light ratio in groups of different size in both 
the bj and rp bands. In constructing the 2PIGG catalogue 
(Eke et al. 2004) and in the present analysis, a well defined 
methodology has been followed based on the extensive use 
of mock catalogues. The starting points are N-body simula- 
tions of the evolution of dark matter in a ACDM universe in 
which model galaxies are added with properties calculated 
according to a semi-analytical model of galaxy formation 
based on the precepts laid out in Cole et al. (2000). The par- 
ent simulation provides a full description of the model galaxy 
distribution, free from the distortions inevitably introduced 
by observational procedures, such as selection effects, obser- 
vational errors, etc. These distortions are modelled carefully 
in order to generate mock catalogues that correspond to 
artificial 2dFGRSs. Since the properties of the parent sim- 
ulation are known, the mock catalogues provide a rigorous 
way to quantify the systematic uncertainties in the derived 
properties of interest as well as a good guide to possible sys- 
tematic errors. Of course, the mock catalogues also enable 
a detailed comparison between the real data and the model 
from which the mocks were generated. 
The main results of this paper are: 

(1) The galaxy luminosity function is different in groups of 
different mass. These luminosity functions are moderately 
well described by Schechter functions for which, as the group 
mass increases, the characteristic luminosity, L„, increases 
and the faint end slope, a, becomes more negative. How- 
ever, the galaxy luminosity functions in the groups found in 
the original simulations are not well described by Schechter 
functions. For example, they exhibit a 'bump' at the bright 
end whose amplitude varies with group mass reflecting the 
relative importance of the central galaxy and the satellites 
in the groups. Hints of bumps are also found in the largest 
2PIGG clusters. 

(2) The median group mass-to- light ratio, T, also varies with 
halo size (which is most robustly characterized by total bj 
halo luminosity). The mock catalogues indicate that, in the 
2PIGG catalogue, T is reliably determined for groups of 
size ranging from that of the Local Group to the richest 
clusters in the survey. Over this range, which extends from 
L bj = 10 10 h~ 2 L to Lj, T = 10 12 hT 2 L , Tj, T increases by a 
factor of 5 whereas T rp increases by a factor of 3.5. At the 
highest luminosities, T becomes roughly constant in both 
bands. The semi-analytical models predict an upturn in T at 
luminosities lower than Lbj ~ 10 10 h~ 2 L at which galaxy 
formation is most efficient. Unfortunately, the 2PIGG cata- 
logue does not contain enough small groups with sufficiently 
accurate estimates of T to locate the theoretically expected 
minimum. 

(3) For the L bj > 3 x 10 11 h~ 2 L objects, the median mass- 
to- light ratio is T(, T = 466 ± 26 (statistical) h M /L , in- 



dependently of cluster size. Assuming that this value re- 
flects the cosmic mean (i.e. that T in these richest clusters 
has converged to the global value) allows a determination of 
the mean cosmic density, Q m . Adopting the bj luminosity 
density inferred from the 2dFGRS galaxy luminosity func- 
tion by Norberg et al. (2002), leads to Q. m = 0.26 ± 0.03. 
(This estimate includes two small corrections derived from 
the mock catalogues as discussed in Section 4.4, and a third 
small correction for incompleteness in the 2dFGRS parent 
catalogue.) This value and its uncertainty are in excellent 
agreement with the values inferred by Spergel et al. (2003) 
from a combination of microwave background data and the 
galaxy power spectrum in the 2dFGRS. 

The agreement of the 2PIGG results with the predic- 
tions of the simulations is impressive, and provides confi- 
dence in the basic picture of galaxy groups tracing viri- 
alized haloes of dark matter. The simulations suggest two 
fruitful avenues for improving upon this work: (i) probe 
smaller groups to test whether the M/L ratio does indeed 
pass through the expected minimum around a group mass 
of 10 12 h' 1 M ; (ii) create catalogues with better resolved 
groups to see if the non-Schechter nature of the theoretically 
predicted luminosity functions can be revealed. Both these 
targets will require new large redshift surveys that probe 
substantially further down the luminosity function than is 
possible with the current dataset. 
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APPENDIX A: CALCULATING GROUP r F 
LUMINOSITIES 

As the 2dFGRS is a 6j selected survey, extra consideration 
is required in order to include the SuperCOSMOS r F band 
data in the analysis. For instance, it is necessary to deter- 
mine an appropriate local flux limit to which the r F band 
data are complete. The term complete is used here with- 
out reference to the well- quantified incompleteness in the 
redshifts measured for objects in the 2dFGRS parent cata- 
logue, and the possible incompleteness in the parent cata- 
logue itself. A local tf band flux limit is determined from 
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b j" r F 

Figure Al. The cumulative distribution of galaxy 6j — rp values 
for all z < 0.12 galaxies in the 2dFGRS (solid line) and mock 
(dotted line) catalogues. 

the distribution of observed galaxy colours. Fig Al shows 
the cumulative distribution of the galaxy bj — tf values for 
the sets of z < 0.12 galaxies in both the mock and the 2dF- 
GRS data. As at least 95 per cent of galaxies are bluer than 
6,j — t'f = 1.5, the local minimum r F band flux limit is set 
equal to r F ,ii m = bj.iim — 1-5. From this figure, one can also 
infer that the bright bj > 14 limit can essentially be trans- 
ferred into the tf band, such that the tf band sample is also 
complete for rF > 14. Thus, a 'complete' rF-limited sample 
can be made to a locally defined flux limit corresponding to 
14 < rF < 6j,iim — 1-5. This lower flux limit removes about 
29 per cent of the galaxies from the sample, whereas only 
~ 0.4 per cent were discarded for being brighter than the 
upper flux limit. 

Having found flux limits defining a complete tf band 
sample, the intention is, as before, to tot up the galaxy lu- 
minosities taking into account the weights associated with 
the redshift incompleteness, and then apply the appropri- 
ate Schechter function correction for the galaxies that are 
too faint to make the flux cut. However, the galaxies in- 
cluded above the rF flux limit are now a bright subsample 
of the 2dFGRS parent catalogue. Thus, they are no longer 
drawn from the same population as the subsample of parent 
catalogue galaxies without measured redshifts, and so equa- 
tion (2.2) is no longer an unbiased way to correct for this 
redshift incompleteness. This complication was sidestepped 
by redistributing only the weights of galaxies without red- 
shifts and with 14 < rF < bj,n m — 1-5. Thus, the subsample 
of galaxies with redshifts and sufficiently high rF band fluxes 
is now a similar population to those galaxies without red- 
shifts whose weights are redistributed. As only about 2/3 of 
the 2dFGRS galaxies satisfy the rF flux limits, the weights 
were redistributed to the nearest 7 projected galaxies only. 
This should match the angular smoothing scale with that 
used in the bj band. 
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Figure A2. The rp band luminosity functions, estimated using 
a 1/Vmax method, for all z < 0.12 galaxies in the mock (dotted 
line) and 2dFGRS (solid line). The dashed line shows the suitably 
normalised Schechter function resulting from an STY fit to the 
2dFGRS data. 



In addition to the correction for redshift incompletion, 
galaxies below the flux limit also need to be accounted for 
in determining the total group luminosity. This was ac- 
complished by assuming that the galaxy luminosity func- 
tion in the groups can be fitted by a Schechter function 
with (M»,a) = (-20.97, -1.16). This is the STY-estimated 
Schechter function for the rp band galaxy luminosity func- 
tion of all z < 0.12 galaxies in the 2dFGRS. In calculating 
this, the following k + e-correction, derived using Bruzual & 
Chariot (1993) models, was used: 



k + e = 



z + 2z z 
l + 9z 2 - 5 ' 



(Al) 



Choosing a normalisation to match the luminosity function 
amplitude inferred using a 1/Vmax estimator leads to = 
(1.2 ± 0.1) x 10 _2 (/i _1 Mpc)~ 3 . This value, which neglects 
the incompleteness in the parent catalogue of the 2dFGRS, 
gives a mean tf luminosity density in the Universe of p rp « 
2.2xlO" 2 /iL Mpc~ 3 . As can be seen in Fig. A2, the rp band 
luminosity functions estimated using the 1/V ma x method are 
similar in the 2dFGRS and the mock catalogue (which yields 
M* = -20.79 and a = -1.17). For the r F band, M = 4.57. 
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